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Abstract. We present a scheme for enhancing the anharmonicity of nanomechanical 
resonators by subjecting them to inhomogenous electrostatic fields. We show that this 
approach enables access to a novel regime of optomechanics, where the nonlinearity 
per quanta of the mechanical motion becomes comparable to the lincwidth of the 
optical cavities employed. In this "resolved nonlinearity regime" transitions between 
phonon Fock states of the mechanical resonator can be selectively addressed. As one 
application we show that our approach would allow to prepare stationary phonon Fock 
states in experimentally realistic devices. Such states are manifestly non-classical as 
they show pronounced negative Wigner functions. We calculate the mechanical steady 
state by tracing out the cavity modes in the weak optomechanical coupling limit and 
corroborate our results by a numerical analysis of the full dynamics including the cavity 
modes. Finally, we show how the negativity of the stationary states' Wigner function 
can be read off the output spectrum of a weak resonant probe laser and analyse possible 
imperfections of the envisioned setup in detail. 



PACS numbers: 85.85.+j,42.50.Dv,42.50.Wk,03.65.Ta 
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1. Introduction 

Substantial progress in fabricating high-Q mechanical resonators with high frequencies, 
as well as recent success in cooling them close to the motional groundstate [H El El E] , 
inaugurate a new research field of manifold fundamental interest [H H E]- The regime of 
very low temperature, where a quantum mechanical description predicts only few quanta 
of mechanical motion, promises potential insight into some fundamental questions of 
decoherence, as well as various technical applications that use the expected quantum 
behaviour [8l |9l [10] . Fundamental questions, concerning the border between the classical 
(macroscopic) and the quantum (microscopic) worlds |TT], trigger a natural interest in 
preparing quantum states of "as large as possible" objects and demonstrating their 
distinct quantum behaviour by appropriate measurements. 

Regarding this major goal, it is important to stress that the dynamics of a 
purely harmonic quantum system is analogous to its classical dynamics, in the sense 
that expectation values of observables follow the classical equations of motion fn\ . 
Therefore, it is a common approach to introduce nonlinearities in a quantum system in 
order to detect quantum behaviour. While there is the possiblity to go to the strong 
optomechanical coupling regime fiSl IT¥[ 115] and make use of the nonlinear nature of the 
standard optomechanical coupling, or to couple to a nonlinear ancilla system [21 [6l [16] , 
we use here a different approach: We introduce a novel kind of optoelectromechanical 
system that features a tunable mechanical nonlinearity. This nonlinearity originates 
from intrinsic elastic properties fTli \TE\ \W\ [20] and is enhanced with the help of 
electrostatic fields. This has the advantage that the cavity sideband driving technique 
is preserved, as this only works in a regime where the optomechanical interaction is 
weak and can be approximated by a linear coupling. The regime of large mechanical 
nonlinearity then enables new means to control the mechanical motion at the quantum 
level, if combined with the coupling to a high-Q optical cavity mode. 

The intrinsic anharmonicity in the mechanical motion of real world micro- and 
nanomechanical resonators is usually very small and therefore only relevant in the regime 
of large oscillation amplitudes. In order to obtain a sufficiently large anharmonicity for 
displacements at the scale of the quantum mechanical zero point motion, we propose to 
use electrostatic gradient forces |21j . These forces originate from the dielectric properties 
of the resonator material inside a inhomogenious, external field. They can be used to 
effectively reduce the resonator's stiffness and therefore its resonance frequencies. In 
turn, this has the effect that the zero point deflection is enhanced up to an extend that 
the nonlinear contribution becomes important. 

Using this technique, the nonlinearity per phonon can be made comparable to 
the linewidth of a high flnesse cavity. This "resolved nonlinearity regime" enables to 
selectively transfer population from phonon Fock state \n) to \n + 1) (|n') to \n' — 1)) 
by driving the spectral line corresponding to n {n') phonons in the blue (red) sideband. 
While the optomechanical technique of driving sideband transitions is well known from 
sideband cooling [221 ES], the anharmonicity of the mechanical motion is a new and 
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crucial feature of our approach. 

In this paper, we exphcitly derive the fundamental mode properties of a nonhnear 
mechanical resonator subject to the aforementioned gradient forces to obtain a suitable 
model for the mechanical degree of freedom. Using that model, we analyze the effect of 
several laser driven cavity modes that are each coupled to the resonator. This enables 
the identification of parameters that lead to distinct nonclassical steady states. 

To demonstrate the capabilities of our scheme, we show results that have been 
calculated considering state of the art experimental components. For the mechanical 
resonator we assume a carbon nanotube (CNT) that is mounted on a chip, together with 
electrodes that generate the electrostatic field. This forms a nanoelectromechanical 
system (NEMS) based on existing devices. For the cavity we consider a high finesse 
toroidal or fibre based microcavity which couples to the CNT via its evanescent field. 

The remainder of the paper is organized as follows. In section [2] we describe the 
nonlinear dynamics of thin rods starting from elasticity theory and derive the resulting 
fundamental mode Hamiltonian. In section |3] we describe our approach to enhance the 
intrinsic mechanical nonlinearity by mode softening with gradient forces. After that we 
introduce the optomechanical model with several laser driven cavity modes in section 
|4j In section [s] a reduced master equation (RME) for the mechanical motion is derived, 
followed by an analysis of its steady state in section [6j We present numerical results 
in section [7] and in section |8] we propose a method to measure the mechanical state 
in an experiment. Finally, we quantitatively discuss central physical properties of our 
setup, namely the optomechanical coupling mechanicsm in section [9} as well as potential 
additional photon losses introduced by the employed electrodes in section [TOj 

2. Elasticity and fundamental mode description for nanoresonators 

The harmonic description of the motion of thin rods is based on only considering the 
bending energy for small deflections |2l]. We consider here thin rods, which means that 
the cross-sectional dimensions, such as the radius for circular cross sections or width 
and depth for rectangular cross sections, are much smaller than the length L. We also 
consider the rod to be homogeneous along the longitudinal axis, here parametrized by 
X G {0,L}, with constant mass line density fi and use thin rod elasticity. The planar 
deflection in the transverse direction is described by a field y{x) where the end points 
at X = and x = L are fixed, i.e. y{0) = y{L) = 0. The Lagrangian within this 
approximation reads 




(1) 



with a kinetic part as well as the bending energy 




(2) 
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R = -W 




u,Ay) = idl' -dl)/dl = -y/R 




Figure 1. Sketch of the deflected rod. a): The local strain u^x is given by the distance 
to the elastic line y as well as the local curvature y" and determines the energy density 
dE/dV = \yu^x- b) Resonator with electrodes modeled by point charges (?,(?'. The 
field profile (red) leads to an inverted parabola for the dielectric potential (blue) around 
the equilibrium position. 



Here, J-" = Y A is the linear modulus of the rod given by the Young modulus Y of 
the material times the cross-section area A. = 4 f iPdA is the ratio of bending 

A J cross ^ ° 

to compressional rigidity and depends on the cross-sectional geometry, where y is the 
in-plane coordinate within the cross-section that is directed along the deflection, with 
origin at the neutral line. For a rectangular cross-section of thickness d, k = d/\/l2, 
whereas for a circular cross-section with radius R, k = R/2. For a cylindrical shell like 
a nanotube, one finds K = R/y/2. The energy l\2n results from the fact that for small 
curvature y", the local strain inside the rod is linear with respect to the distance y to 
the neutral line (see FIG.[T}, while the free energy density is quadratic with respect to 
the strain. The Lagrangian ([T]) leads to the equation of motion 



/ii/ + ^/€V'^=0, 



(3) 



where y^'^'^ denotes the 4-th derivative with respect to x. As this equation is linear in 
y and its derivatives, it leads to harmonic dynamics. The implementation introduced 
later is best decribed by clamped boundary conditions y'{0) = y'{L) = for the rod. In 
that case (|3| has the eigenmodes 



1 

Cn 



sm^Unx/L) — sinh(^'„a;/L) cos{unX / L) — cosh(z/„x/L) 



sin(z/„) - sinh(i^„) 



cos(z/„) — cosh(z/„ 



(4) 



with frequencies a;„ = CgH {un/ L)'^ , where Cg = \fTj[i is the phase speed of 
compressional phonons along the rod. The z/„ are the roots of the transcendental 
equation cos (z/„) cosh (i/„) = 1, i.e. v\ = 4.73. The Cn are normalization constants 
chosen such that max {0„(x)} = 1. We choose this normalization so that the coefficients 
in a mode expansion repesent the maximum amplitudes of the deflection associated to 
each mode. Introducing now the canonical momentum Il{x,t) = 6C/6y{x,t), as well as 
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the expansion of the field into the modes 

y{x,t) = Y,Mx)Mt), (5) 

n 

yields the Hamilton function of a harmonic oscillator for each mode 

with the defiection Xn and mode momentum Vn = m^dtXn for the n-th mode, as 
well as the effective mode masses m* = (j)'^(x)dx. Corrections to this harmonic 
description that lead to nonlinearities originate from a streching effect that occurs due 
to the deflection if the end points of the rod are fixed in space [25] . The resulting strain 
leads to an additional energy 

where the stretched length is = J a/1 + {y')'^dx ~ L + | J dx(y')^ with L being the 
zero deflection length. Including this streching energy leads to a nonlinear extension of 
the Hamiltonian for the fundamental mode 



where we introduced the fundamental frequency Uj^fl and the effective mass of the 
fundamental mode m* ^ 0.3965 fiL. The anharmonicity is given by 

2 



f^(P'\xL)dx 
uf {m* /m) 



The nonlinearity also induces a weak coupling between modes that can be neglected 
for low enough temperatures. We quantize the Hamiltonian ([s]) by introducing bosonic 
mode operators b\ b given by 

(10) 




and obtain the phononic description 



H = huj^,ob^b + h^{b^ + b)\ (11) 



with the nonlinearity parameter Aq = I^^zpm/^) where xzpm = ■\/^/'2'>n*^m,o is the 
zero point motion amplitude of the mode. Here, the frequencies Wm.o and Aq refer 
to fundamental mode properties that result only from the intrinsic elastic forces and 
without applying any external forces to the rod. As we will describe in the next section, 
external forces can be used to tune the resonance frequency — )■ Wm, which will in 
turn change xzpm and Aq — )■ A. 
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3. Enhanced nonlinearity by mode softening 

As stated above, the nonhnearity per phonon scales as A . Therefore, we 

intend to lower the harmonic frequency cUm.o by applying external forces. One way to 
change the mode frequency that has been discussed previously, is to add an additional 
external force along the rod axis, that causes compressive or tensile strain [251 EE] . An 
alternative approach, that promises better control but yields the same potential for the 
fundamental mode, is to use a static electric field that is strongly inhomogenious in 
the direction of defiection. If the rod material shows suitable dielectric properties, this 
produces an additional, inverted square potential with respect to the defiection, see FIG. 

To describe our approach we consider tip electrodes with suitable potentials at 
each side of the nanoresonator that generate a static field. In our calculation we model 
the electrodes by point charges q, q'. The electrostatic energy associated to the dielectric 
per unit length along the rod can be described by 

W{x,y) = + , (12) 

where x, y are the co-ordinates along the resonator axis and the direction of its defiection. 
i^ll X are external field components parallel and perpendicular to the beam axis and ay ^ 
the respective screened polarizabilities. We can expand W{x, y) to second order in the 
displacement y and get an additional contribution to the Hamiltonian that reads 

^ 1 r d'^W{x,y) 



V _T/(i) , T/(2) _ f dW{x,y) 



dy 



y {x)dx , (13) 



where we dropped the displacement independent constant W{x,0) which is irrelevant 
for the dynamics. Inserting the modes ^ we get 

n Ik 

with Fn = J dyW\y^Q(f)ndx and Wik = J dyW\^_^(pi(pkdx. The electrostatic forces Fn 
cause a static defiection of the rod and shift it's equilibrium position. However if the rod 
interacts with the photon fields of a nearby cavity, those fields will also cause a defiecting 
force. For convenience we choose parameters such that these two forces compensate and 
the equilibrium position remains unshifted (see section |4]). 

In general, the contributions of Wik will cause weak interactions between the modes, 
as Wik is not a diagonal matrix. By diagonalizing the Hamiltonian of the rod in the 
presence of electrostatic fields, one can find new mormal modes. For the parameters 
considered in this work, the induced corrections to the mode shapes are however found 
to be negligibly small. We thus end up with the a mechanical Hamiltonian for the 
fundamental mode of the form, 

H^ = ^ + \nfujlX^ + ^X' + F,X , (15) 
2m* 2 4 
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with a reduced frequency ~ g ~ I ^oo In a phononic description this reads 

= hw^b^b + h^{b^ + bf + FoXzpm(&^ + h) , (16) 
where the nonhnearity per phonon A is now increased by a factor C,"^ = (00^^/00^)'^ . 



4. Optomechanical Model 



We consider a typical optomechanical setup with a micro toroid cavity coupled to the 
nanomechanical resonator, where the displacement of the latter modifies the frequencies 
of the cavity modes. The cavity modes with frequencies Ui that contribute to the 
dynamics are described by creation and annihilation operators a^^ and Oj. They are 
each driven by a classical laser of input power Pj and frequency cul.i- The coupling 
strength to the different cavity modes is given by GqjXzpm, where Gqj = ^ is the 
optical frequency shift per deflection deflection X. Thus, in a frame rotating with the 
laser field modes, the Hamiltonian describing the coupled system reads {h ^ 1) 



H 



a'j + aj 



GojXzpMaj'^aj [b^ + b) 



Hrr 



(17) 



'j and Qj/2 



Here we have introduced the laser detunings denoted by = ui^j 
^y PjKex,j/}TiOL,j, with Kex,j being the cavity decay rates into the associated outgoing 
electromagnetic modes. 

Both the field inside the cavity and the mechanical motion are subject to damping, 
which in the regime of weak optomechanical coupling, small nonhnearity and low 
mechanical occupation is well described by a master equation with Lindblad form 
damping terms. With the decay rate rates for the cavity modes nj and the mechanical 
damping rate 7m, this master equation reads 



P 



:i8) 



Here we introduced the Lindblad form dissipator T>{6)p 
as the Bose occupation number n 



2dpd^ — d^dp — pd^d, as well 

exp (^f^ j — 1 of the phonon bath mode with 
frequency at tempertature T. A more precise treatment of the mechanical damping 
accounts for the fact that, due to the mechanical nonhnearity, there is more than one 
bath mode resonantly coupling to the mechanical mode. However, for high-Q resonators, 
(18) proofs to be sufficiently accurate. 

As usual, we expand the cavity field operators around their steady state values and 
adpot a shifted representation aj — )■ aj + aj, with aj 



Qj/{2Aj + iKj), in which the 
master equation has the same form as in ( 18 ) but with the shifted system Hamiltonian 



aj + H.c. ) {b^ + b) 



+ u^b^b+^ib^ + b)^ (19) 



Nonlinear nanomechanical resonators for quantum optoelectromechanics 



8 



with g^j = 2ajGQjXzpM and where we dropped the nonhnear term oc aj^aj(6^ + b) in 
the couphng, which is vahd for {aj^aj) ^ We also assumed that the static electric 
fields are chosen such that Fq = —h GojIctjP, so that the equilibrium position of the 
mechanical oscillator is undeflected. 



5. Reduced Master Equation 



As we are interested in the dynamics of the mechanical resonator, we now choose to 
trace out the cavity degrees of freedom to obtain a reduced master equation for the 
mechanical motion. To do so, we focus on the parameter regime suitable for back action 
cooling. Indeed, we will find that the processes leading to nonclassical steady states 
can be interpretated as a combination of cooling and heating processes with respect to 
different number state transitions. Successful cooling (or heating) requires g^j <C nj and 
^7m ^ 5'mj- The first condition ensures that the cavity modes stay in the vacuum state 
with respect to the shifted represention. This enables us to trace out the cavity modes 
safely. The second condition ensures that the influence of the cavity on the mechanical 
motion is greater than the influence of the thermal environment. We use this separation 
of frequency scales to adiabatically eliminate the photon degrees of freedom. An outline 
of this approach is given in appendix refappendix. 

For our analysis, we focus on the regime of small nonlinearity A <^ and apply 
a rotating wave approximation (RWA) for the mechanical part Hja H[^. This 
diagonalizes the mechanical Hamiltonian in the number state basis (again assuming 



(20) 



with the shifted frequency cu^ = cUm + A' and A' = 6A. Note that in this approximation 
the energy eigenstates of mechanical motion are Fock number states. It is useful to 
switch to an interaction picture with respect to the nonlinear mechanical mode which 
requires to write the mechanical mode operator b in its matrix representation in the Fock 
basis, bn = \fn\n — as those projectors rotate at frequencies that now depend on 

the Fock number n. 



with Sn — uj'^ + X'ln — 1). This Fock number resolved dynamics consequently appears in 
the resulting reduced master equation for the mechanical density operator /i — Tr^ {p} 



-1 



A' 



<b^b + -b^b^bb + J2 iKPnbn + A-^;^bX) , 



n,3 



+ 



(22) 



4E 



V{bn)^i 



J2 + 7mri 
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Figure 2. Level diagram of tlie coupled cavity-resonator-system. \n^n'j) denotes a 
state with n phonons and n' photons in the cavity mode j while the other cavity 
modes are in the vacuum state. |n,„Oc) means all cavity modes are in the groundstate 
(empty). Due to the nonlinearity, the mechanical energy level spacings depend on the 
fock number n and can be addressed individally. Here we have sketched the driving 
of the phonon transition |0) |1) by a blue detuned laser Ai — ojl,! — uji = 5i as 
well as those of transitions |2) |1) and |3) |2) by red detuned lasers A2 = 
A3 = -S3. 



where we have introduced small cavity induced frequency shifts for the nanoresonator 



as well as the transition rates j . The latter are exact analogues to the cooling and 
heating rates recovered in [27] and are here given by 



^L-^?m,,4(^^.^£)2 + ^2- (24) 

In the context of cooling, the analog rates determine the effective temperature of the cav- 
ity modes as a cool bath and so also the final temperature of the mechanical resonator. 
Here, in contrast, the influence of the A^j is focussed on driving transitions between 
mechanical Fock states \n) and |n — 1), where the direction of the drive is given by the + 
or — sign. The magnitude of j as a function of the detuning Aj is a Lorentzian with 
resonance ±(5„ and width Kj. This enables to selectively address specific transitions by 
choosing the detuning of one of the lasers to match the corresponding energy spacing. 
In figure [2| a sketch of this mechanism is given using the relevant energy level diagram. 
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6. Steady state solution 

We assume that the density matrix /i will be diagonal in steady states and seek for a 



solution of ( 22 ) with the ansatz 

^i = J2Pn\n){n\. (25) 

n 

Since within this ansatz the unitary part vanishes, this yields the steady state relations 

' -F{n). (26) 



ra-l 



These relations together with the normalization Yin = ^ define the steady state 
of the mechanical oscillator. It can be calculated numerically by imposing a cutoff 
number uq, for which Yln=no+i ^ ^' ^/^o = T^n Il'i=o -^(^ + 1) ~ 



= 710 + 1 



nr=o -^(^ + -'-) ^^"^ other P„ successively with the recursion relation (26). Note 
that without laser drives the j vanish and we obtain relations characterising a ther- 
mal state at the bath temperature T. If instead we let A — )■ 0, the ratio P„/ Pn-i becomes 
n- in dependent, implying a thermal state at a different temperature depending on the 
laser detunings and amplitudes. 



7. Setup and results 

For the experimental realization we envisage a setup as shown in figure |3j comprising a 
NEMS chip and a high finesse toroidal micro cavity [291 EO] ■ The NEMS chip consists 
of a single-walled CNT (cs = 21000 ms~^) acting as the nanoresonator [311 [32], ^ well 
as electrodes that generate the softening field. The nanoresonator interacts with the 
evanescent field of the micro cavity via dielectric gradient force (see section [9]). A 
CNT is a favorable implementation, as the intrinsic nonlinearity per phonon can be 
found to scale as Aq oc {K^m) ^, which suggests a small transverse dimension and mass. 
Furthermore, carbon nanotubes show ultra-low dissipation [33], that is even decreased 
in the regime of small amplitudes |34j . 

We calculated the steady state of mechanical motion for several sets of laser 



detunings A,-, both by use of the relations (26) and by solving the full master equation 



(18) with the shifted system Hamiltonian ( 



9). In order to obtain self consistency with 



respect to the applyied RWA (20) on the analytical side, as well as justification for 
a Fock number cutoff in the numerical treatment, we aim to produce steady states 
with small phonon occupations ~ 1. We also want to produce states with a 

significant negativity in the Wigner phase-space distribution, as this is a commonly 
accepted indicator of a nonclassical state. Therefore, the most promising choice is to 
aim at a single phonon Fock state |1) [35]. This state's Wigner distribution shows 
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Figure 3. a): Micro toroid cavity with nearby CNT (not to scale), b): NEMS-chip 
witli CNT-resonator and electrodes, c) Schematic top view of the NEMS chips with 
the CNT oscillating in the z-y-plane. d) : Schematic cross-sectional view of the NEMS- 
Chip and the cavity rim (yellow). The NEMS-chip is mounted in the evanescent cavity 
field at a distance d to the cavity rim surface. The nanoresonator is displaced from the 
closest point to allow for an optomechanical coupling that is linear in the resonator's 
deflection. 



a negativ peak of depth Vr(0,0) = —2/71 at the origin of the phase-space, which is 
the maximum negativity attainable for physical quantum states. In order to drive the 
nanoresonator to this state, the plausible strategy is to drive the Fock state transitions 
|0) — )■ |2) — > |1) as well as |3) — i- |2) to obtain net cooling. This is done by choosing 
laser detunings Ai = 5i and A2,3 = —^2,3, repectively. In figure ^ we show results for 
this scenario, comprising to different sets of parameters. 

The system parameters considered in figure |4^), b) can for example be realized with 
a (10,0)-CNT of length L = 1 |am and radius R = 0.39 nm, that is softened by a factor 
C = 4.0. The polarisability of the CNT is a\\ = 143 (AneoA^^, a± = 10.9 (^47reoA^j 
[36] and for an electrode tip distance of 20 nm from the resonators axis, we find 
for the maximum fields at the axis E\\ ^ 1.20-10^ V/m and Ej_ ^ 1.78-10^ V/m. 
The cavity is a silica {ric = 1.44) microtoroid with wavelength 27Tc/uj ~ 1.1 (xm, 
circumference = 1.3 mm and finesse J-'c = 3 • 10^ [37]. With an optical mode volume 
Knodc = 4.0 • 10~^^ m^ and an decay length of the evanescent field /cj^ ~ 0.17 |xm, the 
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0123456789 10 11 n 



cj p 

1 



0.75 
0.5 
0.25 




d) 



H FME 
H RME 

F(m,|1))=0.55 



1 2 3 4 5 6 7 



-0.15 




9 10 11 n 



Figure 4. Steady state results for the given effective paramters, comparing results 
for the reduced master equation (RME) and full master equation (FME). The given 
fidelities F and the Wigner representation correspond to the FME results. Plots 
a),b): a;,„/27r = 5.23 MHz, A727r = 104 kHz, .gj/27r = 20.9 kHz, KjI^-k = 52.3kHz. 
The ambient temperature is T = 20 mK and Qm — 5-10®. Plots c),d): oj^n/'^T^ = 
2.13 MHz, \' /2tt = 42.7 kHz, gj/2n = 12.8 kHz, Kj/2tt = 64.2kHz, T = 30 mK and 
= 1.5-106. 



given coupling strength can be reached by placing the NEMS chip at distance d = 50 nm 
from the cavity rim and with a launched input power Pin,j = 0.14 W for each laser 
(assuming k^^/k = 0.1). 

In turn for example |4]:), d) we have L = 1.7 |^m and ( = 3.3 which corresponds 
to E\\ ^ 5.3 ■ 10^ V/m and E± ^ 1.0 ■ 10^ V/m, as well as = 1.80 mm, 7^ = 2- 10^ 
Knode = 5.5 ■ 10^^^ m^ and = 2.6 mW. Here, we also assumed a decrease of the 
mechanical Q by the softening about a factor C to Q = 1-5 • 10®. 

Of course, the presented scheme also allows for the preparation of other Fock states, 
or mixtures of Fock states. A few examples are given in figure [5] 

8. Measurement via output power spectrum 

The steady state of the nanoresonator can be probed with an additional laser, weakly 
driving one cavity mode on resonance, i.e. A = 0. This mode has to be distinct from 
the modes used for the state preparation. Then, the occupation of different Fock states 
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H FME 
H RME 

|2)) = 0.86 - 
Ai = 62 
A2 = -^3 
A3 = -^4 



1 2 3 4 5 6 7 



9 10 11 n 





H FME 
H RME 

F(/i,|3)) = 0.80 - 
Ai = 63 
A2 = -54, 
A3 = -^5 



1 2 3 4 5 6 7 



9 10 11 n 



1 2 3 4 56789 1011 n 

Figure 5. Steady states for a nanoresontor with the parameters given in figure |4 
b) and different sets of laser detunings Aj . 



can be extracted from the sideband structure in the power spectrum 

00 

S{u) = ^ [ dre-^('^-'^L)^(aoutHi + ^)«out(r))ss. 
Ztt I 



(27) 



The output field aoutit), is related to the intra cavity field a(t) via the standard input- 
output relation [3S], ^out = y/i^o, + ajn, while the dynamics of the intra cavity field can 
be described by a quantum Langevin equation 



d = --a - iGoa^zPM (^"^ + b) {a + a) + A/K^5ain(i) + - f^cxScinit) , (28) 

where 6ain, Scin are the fiuctuations of the input field in the laser mode and the other 
bath modes. Defining a^'^^ to be a solution for Go = 0, we can integrate (28) formally 



and apply a Dyson series type expansion to first order in the coupling strength per 
photon, Go^^zPM, to find for the motion of the cavity modes 

t 

ait) ^ a(°)(t) - iGoXzPM J dre-«/2(*-r) [^t(^) + [^(o)(^) ^ 





bnjt) 

i6n + k/2 



+ H.C. 



(29) 
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1 2 3 4 5 6 7 8 9 1011 -^t -^s -S2 Si 5i 62 S3 6i 



Figure 6. Left: Example result. Right: Corresponding output spectrum around 
the upper and lower sideband of a resonant driving laser. The peaks correspond to 
different diagonal entries of the mechanical density matrix. They are separated by A'. 
The relative weights of Stokes and anti-Stokes peaks reflect the ratios P„/P„_i. 



Here, the contribution of the fluctuations is included in the free fleld solution a^°^(t) 
and the input fleld operators have already been written in a shifted represention 
o-inif) — 5ctin(^) + (c^in) 55 which splits off the coherent part of the input. 

The spectrum consists of a contribution at the laserline u = wl, as well as a 
sideband structure due to the interaction with the mechanical mode. We concentrate 
on the contributions to the sidebands, which to lowest order in gm/'^m are given by 
the flrst order two-time correlations of the mechanical motion. The latter can be 



calculated with help of the reduced master equation (22) and by use of the quantum 
regression theorem. The only nonvanishing contributions in the steady state are 
{biit + r)6„(t))ss = eK-^"«/2]-nP„ and (6„(t + r)6t (t))ss = el-''--^"^/'VnP^_„ where 
we neglected the contributions of the A^"^ as they induce only a small shift of the peak 
positions. Using {A{t — T)B{t)) = {A{t + r)fi(t))* we flnally obtain 

(30) 



I 1 2 

S[Uy^UL) = — y — — — 



%ff r^, . 1^ I %fT 



with 6u = w — cjl, the distance from the input frequency. Hence the spectrum 
shows Lorentzian peaks, separated by frequencies 5n that match the mechanical energy 
level spacings from the laser line. The linewidth of the peaks is given by 7"^ = 
n{A'l + Al + -f^[2n + l]} + {n-l){A''_-'+im[n + l]} + {n + l) {Al+^ + -f^n} , which 
satisfles 7"^ ^ A for small Fock numbers n provided that n7m ^ Qm- Making use of 
expression ( pO) ), the ratios Pn/Pn-i can be measured by comparing integrated peak 
intensities /(w) at positions w — = ±5n as 

/(wL-^n) Pn-l' ^ ^ 

These ratios need to be measured up to a Fock number above which the spectral lines 
become invisible indicating that the occupations probabilities P„ become vanishingly 
small. Here, the stationary nature of the motional states considerably simplifles their 
measurement as compared to phonon Fock states that might be prepared by projective 
measurements SOI SI] • 
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In the following two sections we analyze important practical aspects of the 
envisioned implementation with a carbon nanotube coupled to a toroidal microcavity in 
more detail. Readers who are only interested in the results of the discussed mechanism 
may thus directly turn to the conclusions section. 



9. Optomechanical coupling via the evanescent field 



In this section we derive an estimate for the coupling between the mechanical 
displacement and the cavity field. The coupling arises as the energy of the dielectric 
oscillator in the evanescent electric field depends on the field strength at the location of 
the oscillator. As the evanescent field decays with distance to the cavity rim, altering 
the oscillator-cavity distance by displacing the oscillator results in a change of energy. 
Thus, the coupling part of the Hamiltonian is given by the dielectric energy 



P{r) ■ E{r) 



(32) 



where the polarization is P (f) = V ■ E{r) with "V?" the screened polarizability tensor 
and E the external field. The integration is taken over the oscillator volume. 

Given the dimensions of the microcavity, the cavity rim can be locally modeled 
by a straight, cylindrical fiber. We introduce cylindrical coordinates (r, 0, 2;) with the 
5-direction along the fiber axis (compare FIG. ^) and consider TEo,n modes, as a trans- 
verse electric field is advantageous with repect to possible losses that are discussed in 



section 10 The corresponding transverse fields are given by 



r < ar 



r > ar 



Br 



— 1 



.fc|| dB, 

j2 Qj. 



Eu, = — -z—Br 

k\\ 



(33) 
(34) 



Br 



— 1 



h dB, 
dr ' 



Ein = — -;—Br 

ku 



(35) 
(36) 



where axial field reads 

B,{r, z) = B-MUiry^^^' , r<a,, (37) 

i?,(r, z) = ^BM%j^e'''' ,r>a,, (38) 
Kq '{K±ac) 

with the modified Bessel function K^^^ and the Bessel function of first kind Jo, and 
where k± and 7 are the transverse wave vectors outside and inside the waveguide, 
respectively. The ratio of the field maximum inside the waveguide to the field at r = ac 
is given by ^ = Bz{ac)/ Bz{0) ~ 0.4 if ac7 coincides with the first zero of Ji{x), which 
is ~ 3.8. We further approximate the evanescent electric field for r > Oc, using 
2K±ac ^ 1 and assume a frequency well above cutoff, so that k\\ ~ Uck ^ l/oc. With 
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these approximations, the evanescent field can be written 

EJr) ^ -i— eX(0)./^e-"^(^-"^)e'"'=^ . (39) 
K± y r 

This field will later be used to estimate losses induced by the electrodes. In order 
to determine the optomechanical coupling strength, we write the electric field in its 
quantized form 



= \r^i^^ + ^Wi^)^ where U{r) ^ ]lL \ —^e~ 



K,j_{r-ac) 



(40) 



with photon operators a, a''' and the mode function U{r) for the evanescent part. Here, 
Knodc ~ 0.5 TcalLc is the cavity's optical mode volume and ^ denotes the ratio of the 
maximum field inside the waveguide to the field at the waveguide surface. The zero 
point motion xzpm of the oscillator is small compared to the evanescent decay length 
1/k,± of the evanescent field and the same holds for the transverse dimensions of the 
oscillator, so that the electric field is assumed to be constant everywhere inside the 
integration volume. Thus, by linearizing E{r) about the equilibrium position of the 



nanoresonator and comparing the Hamiltonians (32) and (17) using (40), we find for 
the opto-mechanical coupling rate, 

Go^3.o "-"f"f""^' e-^-^^Ceo.. (41) 

Here, we neglected the contribution of the perpendicular polarizability since, for carbon 
nanotubes, the perpendicular polarizability is typically one order of magnitude smaller 
than the parallel one. The geometry of the setup we consider here is illustrated in figure 
[3] and the dependence of alignment and positioning of the nanotube is accounted for 
in the correction factor Ccorr- For a TEo,n-iiiode of the cavity field, the electric field 
is directed along ip, i.e. tangential to the cavity rim, while the nanotube is aligned 
along Cnt = cosO'z + sin 9'y. In addition, the defiection of the tube is not aligned 
with the electric field gradient. The direction of the defiection and the direction of the 
field gradient are given by Cosc = cos 6'y — sin 6'z and f = cos ipx + sin ipy. Finally, if 
d is the distance of the chip to the cavity rim, the actual distance of the oscillator 
to the rim is d' = (oc + d)/ cos 99 — Oc. Taking together all these corrections, we 
find for the correction factor Ccorr ~ e"^''^'^'^"'""'^)^'''^'^'^"^) sin^ ^' cos6'' cos^ (^sin(^. This 



correction factor is maximized for sin^^^ = 2/3 and (p^ ^ 1/ a/2k_i_(c? + ac), resulting in 
1/Ccorr ~ 15.5 for a refracting index Uc = 1.44. 



10. Cavity losses induced by metallic nanotube electrodes 

Preparing a nonclassical steady state demands a high-Q optical cavity as the linewidth 
K needs to be comparable to the mechanical nonlinearity A. Using conventional metallic 
electrodes for the mechanical frequency softening can potentially increase the cavity 
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losses. To achieve the necessary low losses it is crucial to have deep subwavelength 
transverse dimensions for the electrodes. 

We now give an estimate of the photon losses that may be induced by such electrodes 
and show that these are negligible. To do so, we model the electrodes as metallic 
cylinders and assume that their radius is much smaller than the wavelength of the 
evanescent cavity field. We assume the electrodes to be parallel to the waveguide, with 
a small misalignment angle 6. For the considered TEo,n modes, the losses result from the 
misalignment only, as for 6 = there is no field along the electrode axis. The resulting 
finesse J-" for the cavity is composed of the finesses J-i associated to the different decay 
channels, which by assuming independent loss channels reads 



For each loss channel, the finesse J-^ is determined by the ratio of the incident power 
and the fraction of power that it lost through the respective channel. The loss channels 
that will be considered subsequently are (1) scattering by the "bulk" of the electrodes 
modelled as a single metallic cylinder , (2) scattering by the "gap" between the electrodes 
and (3) absorbtion. 

10.1. Incident field and incident power 

For our calculations, we introduce new cylindrical co-ordinates {r',(p',z') for the 
electrode with the z-direction along the cylinder axis. We first express the field that is 
incident on the electrode in the electrode coordinates. The result will later be used to 
determine the scattered and absorbed fractions of the incident power. 

To express the incident field on the electrodes in the electrode coordinates, we 
consider the projection of the electric field vector on the electrode axis (see FIG. [7]) 
and assume that the field is constant over the electrode cross section with dimension 
i? -C kJ^. This yields 



We only need to determine the parallel field component E^,{z') here, as the 
other components in the relevant scattering field modes are then already completely 
determined by the wave equations for the field at the electrode. The origin of the 
primed axis lies at {—d — a^, 0, 0) and the directions translate as 




z ■ ipE^ 



(43) 



z = sin 6y + cos 6z 
= — sin {px + cos ipy . 



(44) 
(45) 




(46) 



(47) 
(48) 



z 



z' COS 6 . 
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Figure 7. Relative alignment of electrode coordinates (r, (p, z) and waveguide 
coordinates (r', (^', z'). The misalignment is measured by a small angle 9. 



For our further calculations it is convenient to express the incident field via its Fourier 



transform Ef^\k') = J^^ El'^\z')e-'''''' dz' . Using this and equations (|44|)-(|48|) in Q 



.(in) 



and (39) we can express the incident field as 



^a,{d + a,)^e-''-'BM I dx ^^^^,^3^, , (49) 



F 



where we have substituted x = z'sm6/{d + Oc) and k = {k' — nckcos6)/{K±sm6). 
One can find an approximation for the integral F by applying the method of 



steepest descents, using K±{d + Oc) ^ 1 and d ^ k]_ <^ ac, which yields |F| ^ 
I {d + a(^K\_Q~'^^^'^^°"=^^^^ for l^^l <^ 1 and is maximized by k' = k for small angles. 
Finally, we want to estimate the incident power, which is given by 



dA-\S\ 



27rncC 



rJi('yr)dr ^ 0.51 



CUr 



ka^ 

7 



(50) 

where we have used the Poynting vector S = E^Br/]iQ with the vacuum permeability 
]1q. Now that the incident fields are formulated, we continue to calculate scattering and 
absorption losses. 



10.2. Scattering losses 

We estimate an upper bound for the scattering losses by modelling the electrodes 
as perfectly conducting cylinders, as this maximizes the scattering. We expand the 
scattered field into solutions of the wave equation in cylindrical coordinates for the 
electrode. As the radius of the electrode R is much smaller than the wavelength of 
the incident field A, all contributions are supressed like (R/X)^, as can be seen by a 
Taylor expansion of the corresponding cylindrical harmonics, except s-wave scattering 
of TM modes, for which the supression is only logarithmic. We further assume that the 
incident field is constant for a given cross section and determined by the field at the 
electrode's center. Moreover, we neglect multiple scattering between the electrode and 
the waveguide. 
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We determine the scattered field by use of tlie boundary condition Ez'{r' = 
R, ip', z') = + i^^f-* = at tlie surface of tfie electrode. The transverse fields of 
TM solutions for the scattered field with 2;'-dependence exp {ik'z') are given by 

E'l' = i^^^V.E?) , (51) 

H^I^ = ceo^z' X , (52) 

where k is the wave vector in vacuum and El^\k') = HQ^\K±r')e^'''^' , with = k^ — k'"^ 
and k''^ < k"^. For /c'^ > k"^ the solution is evanescent and does not contribute the 
scattered power. Thus the scattered field can be written 

k 

^^'^"'^'"^ 7 2vr^^' ^^^^ln(v/F^i?)' ' ^^^^ 

where E^^f'{k') is the Fourier transform of the incident field E^^^\k') and where we have 
used the approximation Hq^\x) ~ i^ Inx for x <^ 1. We calculate the scattered power by 
integrating the scattered energy flux across a cylinder with infinite radius -R* — > oo that 

enwraps the electrode, where we use (f' ■ S) = l^{f' ■ {E x H)} = l^ij^^^E^^f}, 

as well as Hq^\x) ~ ^/2/ttx exp [i(x — 7r/4)] for x ^ 1. Thus the scattered power reads 



k 

Ps ~ -ceo / dk' 



4 J (P-fc'2)ln2(VA;2-fc'2)i? 

— k 

k 

<7Tceomaxl E^r\k'] 



2 / dk'- ^ 



\k'\<k 



(A;2 - k'^) In' [(fc2 - fc'2)i?2] 







G 

The Integral G can be estimated for kR <^ 1 which yields G ~ 1/(2 1 ln{2kR)\). So by 



comparing this to the incident power (50) and using sin^ ^ 9, cos 9 ~ 1 for \9\ <^ 1, we 



find an lower bound for the finesse associated to scattering losses 



^■> '"''^ l..(A)e'-^- (54) 



7r^2 \4:TtR 

Here we have introduced the denhancement factor ^ = and assumed d = 1/{2k,±) 
and Oc = A. These assumptions are justified since for typical parameters one has d <^ 
and ac ~ A. For A = 1.55 |J.m, R = 10 nm, Uc = 1.44 and 9 < 10° we thus find 
J-'s > 7 X 10^^. Hence losses due to scattering are cleary sufficiently suppressed for 
electrodes of sub wave length radius that are aligned with 9 < 10° with the cavity rim. 

10.3. Gap contribution to scattering losses 

So far we have considered scattering from one single continuous electrode. Actually, 
our setup comprises two electrodes (see FIG. [s]), separated by a gap in which the 
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nanomechanical oscillator is positioned. We now estimate an upper bound to the 
additional losses that this gap may induce. To do so we model the gap of size D 
as a radiating dipole that is induced by the incident field. We approximate this dipole 
by a perfectly conducting sphere with radius D/2 to optain an upper bound for the 
induced diploic. Thus we find p ~ |eo-D^_E'^"'' cos^^ [12] and for the scattered power, 

/ (■ \ 2 



P, < -eock^D' (^^ 15.(0)1 ) e-^--'^ cos^ . (55) 



Along the same lines as before, using (50), the associated finesse reads 



and shows that scattering losses due to the gap between both electrodes are negligible. 
10.4- Absorption losses 

We consider transparent electrodes with 2D conductivity a that absorb the power 

oo 

P. = ll dz'3?{r4^)} , (57) 

— oo 



since -E^™^ induces a current / = 2tx RaE^^f' in the electrode. By use of (43) and (50) 
the absorbed power reads 



oo 



,2 r „-2K_L(rf+ac)(v/l+S2-l) 



7 y (1 + x'^y^i'^ 



J 



where we have again substituted x = z' sm9/{d + a^) and J ~ {^±{(1 + ac)/n)^^^'^ for 
2K±{d + ac) ^ 1. Thus using d ac yields for the finesse associated to absorption 
A. ~ TTo-flg I _!^^-2K^d g-^ Q Nanotubes exhibit an upper bound in the conductivity 

o'max = 8e^/h pS]. By assuming an off-resonant o = aamax with a < 1, we get 

WnaFe^y 71 RsinO' ^^^^ 

with the fine structure constant ap ~ 1/137. Plugging in the same values as before 
except that now < 3° and R = 2.5 nm and assuming a < 1/20 we find ^ 7 - 10''. 
Hence even though absorption losses clearly dominate over scattering losses, their effect 
can still be neglected for electrode radii R < 3 nm an alignment angles ^ < 3°. 
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11. Conclusions and Outlook 



We have introduced a scheme to novel regime of optomechanics where the 

motion of the nanomechanical resonator becomes anharmonic and thus allows to explore 
genuine quantum dynamics. In our approach inhomogeneous electrostatic fields are 
applied to the nanomechanical resonator to enhance its anharmonicity per phonon 
until it becomes comparable to the linewidth of a high finesse optical cavity. An 
optomechanical coupling to driven modes of this high finesse cavity allows to selectively 
address transitions between phonon number states of the mechanical resonator and 
hence to control the mechanical motion at a Fock state resolved level. The scheme set 
forth thus opens new avenues for optomechanical manipulations of nanobeams, such as 
the preparation of phonon Fock states or optomechanically assisted approaches to the 
buckhng instabihty pH25l[26]. 
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13. Appendix A 



We consider the Liouvillian associated to the master equation (18) in the shifted 
representation (19) 



c{t) = c'Co + CCiiCt) + C2{Ch) 



(60) 



where we have introduced a formal dimensionless parameter ( that labels the different 
involved frequency scales and will later be used as a expansion parameter. After 
changing to an interaction picture for the mechanical mode (including the nonlinear 
part), the three parts of ([60|) read 



(61) 



which describes the free cavity dynamics. 



(62) 



which describes the photon-phonon interaction, where -^i^jP = [-^j^Lp]; ^in,jP 
[Xjbn, p] and with the cavity quadratures Xj = {a*aj + ajaj'^) /\/2\aj\, as well as 



A(C't)p 



7n 



(n 



+ l)e^"("-™)«'*P_p + ne-'^("-"^)^'*P+^p 



(63) 
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which describes phonon dissipation with damping operators V^^p = 2hnph\n — b\rfinP — 
pbl^bn, T^tmP = '^H^pbm — bmbnP — pbmb\,- We now define projection operators 



Vp = Tre{p} ® p 



(th) 



Q = l-V, 



(64) 



with the trace over all optical degrees of freedom of the cavities Trc{p} 
Trc^i . . . Trcj . . . {p}, and use thermal states for the cavity modes 



P 



(th) 



p^^j^ with p^^"^^ 



oo 



|?T,) {n\ 



(65) 



where n^j = n{uj) is the Bose occupation number for the cavity mode with frequency 
ojj. We will use n^j = in the calculations, which is a excellent approximation for 
optical frequencies at cryogenic tempertures. 

We derive coupled equations of motion for the relevant part Vp and the irrelevant 
part Qp, solve the irrelevant part formally and insert it into the relevant part to obtain 
a closed equation for the latter which reads 



Vp = VC{t)Vp + VC{t) / drT^ 



Jo*dr'S£(r')t 



-/;dr"Q/:(r")Q 



£(r)Pp(r). (66) 



As Vp is a stationary state of Cq it follows that CqV = VCq = 0. Furthermore we have 
VC2{C^t) = C2{C'^t)V. We are looking for an expansion in 1/C for C ^ 1, so we can 
expand the time ordered terms to get 



J^dr'QC{r')Q 



£dT"QC{r")Q 



,eQCoQ{t-T) 



l + O 



(67) 



Using this and substituting r' = ( {t — t), we arrive at 



Vp = V [CMCH) + C^iCH)] Vp+VC,iCH)Q I dr'e^-'Q£i {(H - r') Vp{t-T' /e)+0{l/e 



To obtain now the lowest nontrivial order in 1/C, we apply the limit C, oo. In this 
limit all fast rotating terms drop out which amounts to a RWA. Inserting (62) we get 

Vp = VC^Vp + P£SS / dr'e^^-'e^"-' QC[-2,Vp + H. c. , (6J 



If we now insert (62) and apply the projection (64) we obtain 



oo oo 

VC['-\Q j dr'e(^'^"+^«)^'Q£^Pp = -y" dre'^"^Tr, { [X.fot , e^"- p ® pf )]] } 



(69) 
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Here we have introduced the reduced density operator for the mechanical mode = 
Trc{p}. The operators belonging to the mechanical mode can be pulled out of the 
trace, and so by defining the cavity correlation functions 



oo 

G,{u) = 6,, j dre^"^Tr,{X,e^«^X,pf)} 



(70) 



equation (68) becomes 



A = £2/i + ^^^{G,(5„) [hlhnll\-G*{-5n) [6L/i&n] +H.C.} . (71) 



Note that there are no cross correlations between different cavity modes, which is 
expressed in the delta function 5ij. The cavity two time correlations for each mode 
follow the dynamics of £o and therefore 



G,{uj) 



1 



-2i [u + Aj) + hij ' 



(72) 



Rearranging ( 71 ) and switching back to Schrodinger picture yields the reduced master 



equation (22) with 



2 2 



(73) 
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